Creation of macroscopic superpositions of flow states with Bose-Einstein condensates 
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We present a straightforward scheme for creating macroscopic superpositions of different super- 
fluid flow states of Bose-Einstein condensates trapped in optical lattices. This scheme has the great 
advantage that all the techniques required are achievable with current experiments. Furthermore, 
the relative difficulty of creating cats scales favorably with the size of the cat. This means that this 
scheme may be well-suited to creating superpositions involving large numbers of particles. Such 

\Q , states may have interesting technological applications such as making quantum-limited measure- 

f"^ ■ ments of angular momentum. 
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I. INTRODUCTION 

. Quantum mechanics allows objects to exist in a coherent superposition of different states. Th does not depend on 
the size of the system and so offers the fascinating prospect of being able to create superpositions of macroscopically 
distinct states. Multiparticle superposition states have been observed in a number of systems including three photons 
0, Ceo molecules Q, and the internal state of four 9 Be + ions 0. Experimental signatures of larger scale quantum 
phenomena were shown when Rouse et al. Q observed resonant tunneling between two macroscopically distinct 
states in a superconducting quantum interference device (SQUID). The observed superposition was between states 
I/"") ' of different flux or opposite currents flowing around a loop. These currents consisted of approximately 10 9 Cooper 
. pairs, meaning tunneling between two macroscopically distinct states had been achieved. Similar systems have also 
' been used to create macroscopic superpositions of these two supercurrent states 0, |(| . 

Bose-Einstein condensates (BECs) are a promising system for realizing similar results. They are composed of up 
to 10 9 atoms with a high proportion in the same quantum state and are sufficiently cold to undergo a quantum phase 
Oh transition from a superfluid to a Mott insulator |jj . There have already been a number of theoretical proposals for 
producing cat states with BECs in a range of different set-ups 0, El EI El 0] • 

In this paper, we present a scheme for producing a macroscopic superposition of different superfluid flow states (or 
equivalently angular quasi- momentum states) in a ring of coupled BECs [T^ . This is important because it provides 
a direct manifestation of quantum mechanics at the macroscopic level in a new system. As discussed by Leggett 
|13| . such states are important for testing the limits of validity of quantum mechanics. All the steps in our scheme 
are straightforward and should be achievable with present experiments. The BEC system may also have significant 
advantages over SQUIDs since it is highly controllable: the coupling between condensates and the strength of the 
interactions between atoms can be tuned over many orders of magnitude. This scheme has the added advantage 
that the degree of control over the system that is required to create a cat state scales favorably with the number 
atoms involved. This suggests that such a scheme may be well-suited to creating 'large' cats. Finally, we discuss 
how a macroscopic superposition of different superfluid flows may also be of technological interest. One possibility is 
quantum- limited measurements of angular momentum or, equivalently, ultra-precise gyroscopes [Til ]. 
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II. THE SCHEME 



The system we consider consists of condensed atoms trapped by the optical dipole force in a lattice formed by a 
standing wave of laser light. In particular, we shall consider three lattice sites in a ring geometry (see Figure 1). 
Each of the lattice sites is coupled to its neighbors by quantum mechanical tunneling through the potential barriers 
separating them. 

This system can be described by the Bose-Hubbard Hamiltonian [l^ 

H=-J (a^b + tfe + c* a + h.c.) + - (a t2 a 2 + 6 t2 6 2 + c t2 c 2 ) , (1) 

where a, b, and c are the annihilation operators for an atom at the respective lattice sites. The strength of the tunneling 
between sites, J, can be adjusted in experiments by changing the intensity of the standing wave, thus altering the 
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potential barrier between sites. The interaction strength between atoms, U, is at best only weakly dependent on the 
potential, but can be controlled by using Feshbach resonances. In the scheme presented here, we will be concerned 
only with changing the value of J. 

Using Bloch's theorem, which describes the wave function of a particle in a periodic potential, we can write the 
wave function of a Bosc-Einstcin condensate in polar coordinates as, 



ip(r, 9 + 27m/3) = e i2win/3 ip(r, 0), 



(2) 



where n is an integer and t is some number to be determined. We can interpret £ as being the angular quasi-momentum 
of the mode in units of h. Normally when using Bloch's theorem (e.g. to describe electrons moving in a periodic 
lattice) periodic boundary conditions are put in by hand as an approximation to simplify the analysis. In the geometry 
considered here, periodic boundary conditions are true physical constraints needed to ensure that the wave function is 
single-valued everywhere. This means that the description of this system using Bloch states is exact. The boundary 
condition can be written as, 



ip(r,9 + 2ir) = tp(r,i 



(3) 



and we see from (01 that this requires that £ is an integer. Another way of say this is that the angular momentum 
must be quantized in units of h. Moreover, we see from @ that changing £ by an integer multiple of 3 does not 
change the wave function. This means that a complete basis is formed by restricting £ to the values £ £ {0, ±1}. The 
annihilation operators for this complete quasi-momentum basis are, 



(a + b + c) 
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(a + e * 27r/3 & + e i47r/3 c) 



(4) 
(5) 
(6) 



Let us now consider the details of the cat creation process. The great advantage of the scheme presented here is that 
it is straightforward to carry out and within reach of present experiments. Macroscopic superpositions can be created 
simply by starting with the ground state in the lattice when the barriers are low, rapidly increasing the intensity of 
the optical lattice, and then waiting for a certain time before rapidly decreasing the intensity of the optical lattice 
again. 

For simplicity we begin by considering a system of only three atoms |16) (we will consider larger numbers later in 
the paper). We want to create a state of the form, 



|V) = ^= (|3, 0, 0) Q/37 + |0, 3, 0) Q/37 + |0, 0, 3) Q/37 ) 



(7) 




FIG. 1: Diagram of the set-up. Three BECs respectively denoted a, b and c are trapped in an optical lattice in a ring geometry. 
It is convenient to describe this set-up in polar coordinates (r, 9). 
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FIG. 2: Probability distribution of the number of atoms in lattice sites a and b for the initial ground state of the system, i.e. 
all atoms in the a quasi-momentum mode, for N = 3 ©. 



i.e. a superposition of all the atoms having no angular momentum and all having one unit of angular momentum 
clockwise and all having one unit of angular momentum anticlockwise. 

The Hamiltonian for the condensate in the lattice when the potential barriers are sufficiently low that the energy 
associated with tunneling dominates the energy associated with interactions, J ^> £7, is given approximately by 

H = -J (a^b + tfc + c)a + h.c.) = - J (2a f a - /3 f /3 - 7*7) . (8) 

The eigenstates of this Hamiltonian are the quasi-momentum modes @ - @ an d the ground state is when all the 
atoms are in mode a (i.e. the zero quasi-momentum Bloch state). For three atoms, this can be written in the basis 
of the number of atoms per lattice site as, 

IV.) = J_( a t + & t + c t) 3 | W ) 




FIG. 3: Final number distribution of atoms in the a and /3 modes after the cat creation process when N = 3. 
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= -L (|3,0,0) + |0,3,0) + |0,0,3» + ^|1,1,1) 

+ | (11,2,0) + |2,1,0> + |1,0,2> + |2,0, 1> + |0,1,2) + |0, 2, 1» . (9) 

Figure 2 shows the number distribution of atoms in the lattice sites, P(N a , Nb) = \{N a , Nb, 3 — N a — Nb\4>)\ 2 , for this 
ground state. As we might expect, we see that it is relatively unlikely to find all the atoms in a single lattice site, but 
more likely to find them equally distributed between sites. 

Our claim is that, if we start with a state of the form of |J3J, we can create a cat state of the form of by first 
rapidly increasing the intensity of the optical lattice, then waiting for a certain time, and finally rapidly decreasing 
the light intensity again. We now study the details of this process. 

The first step is to rapidly increase the intensity of the trapping laser light, thereby increasing the height of the 
potential barriers between the sites and so decreasing the coupling, J. We want the rate at which this change is 
made to be much greater than rate of tunneling between lattice sites, but still adiabatic with respect to the spacing 
of energy levels at each site. This separation of timescales is possible and has been experimentally demonstrated 
for condensates in an optical lattice |17j . If the coupling between the sites is decreased sufficiently by the increase 
in intensity that the coupling energy is small compared with the energy due to the interactions between atoms, the 
Hamiltonian can be written approximately as, 

F = |( a tV + 6 t2 6 2 + ctV). (10) 

Immediately after the barriers have been raised, the state is given by I© since it has not had time to evolve. If we 
now allow this state to evolve under the influence of the Hamiltonian i|l(J|) . after time t, the state is, 

hA) 2 = -l= e ^(| 3 ,0,0) + |0,3,0) + |0,0,3)) + ^|l,l,l) 

+ ie 4C/t (|l,2,0) + |2,l,0) + |l,0,2) + |2,0,l> + |0,l,2) + |0,2,l». (11) 

We now want to see whether this evolves to a cat state. For a state of the form of Q, we need the the probability 
for all three atoms to be found in mode a, P Q (3), to be one-third and similarly for the probabilities for modes (3, 
P/j(3) and 7, P 7 (3). Calculating these probabilities directly we get, 

P a (3) = ^[41 + 24cos([/i) + 12cos(2C/i) + 4cos(3t/<)] (12) 

P /3 (3) = P 7 (3) = [41 - 12cos(E/t) - 6cos(2[/i) + 4cos(3£/t)] . (13) 

If we pick the evolution time to be t = 2w/(3U), we get P Q (3) = Pp(3) — P 7 (3) = 1/3. This is precisely what we 
want: an equal superposition of all the atoms in a, and all in /?, and all in 7. A plot of the full number distribution 
in the quasi-momentum basis is shown in Figure 3 and confirms that a cat state of the form of Q has been created. 

Finally we rapidly lower the potential barriers again. As before, this is done much faster than the rate of tunneling 
between the potential wells, but slowly with respect to the spacing between energy levels within each well. This leaves 
us with the state Q evolving due to the Hamiltonian (JSJ. Since this Hamiltonian has no terms that couple 
different quasi-momentum modes, the population in each does not change with time and so the cat-like structure of 
the state is preserved. The relative phases between terms, however, do change with time since the different quasi- 
momentum modes have different energies, i.e. the non-rotating mode, a, has a lower energy than the rotating modes 
(3 and 7. 



III. LARGER NUMBERS 



We have now demonstrated how the cat creation process works for three atoms. In order for this to be a truly 
powerful technique, we need to demonstrate that it also works for larger numbers of particles, A > 3, to create states 
of the form, 

\i>) = -±= (| A, 0, 0) a/37 + |0, A, 0) Q/37 + |0, 0, N) afh ) . (14) 
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For simplicity, let us suppose that the total number of atoms in the system is a multiple of three, i.e. N = 3n, where n 
is an integer. The initial ground state of the condensed atoms trapped in the lattice when the energy of the coupling 
term dominates, J 3> U, is 



N N-p 



1 / ^ ! 



\p,q,N-p-q). (15) 



The probability distribution of the number of atoms in lattice sites a and b is plotted in Figure 4 for TV = 30. 

We now follow exactly the same procedure as for three atoms: rapidly raise the potential barriers, hold for time t = 
2n/(3U), and then rapidly reduce the barriers to their initial height. It is straightforward to calculate the probabilities 




FIG. 5: Final number distribution of atoms in the a and (5 modes after the cat creation process when N = 30. 
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for all the atoms to be in a single quasi-momentum mode and, as before, we get P a (N) = Pp(N) = P-y(N) = 1/3. 
The full number distribution is shown in Figure 5 for N = 30 and confirms that a cat state of the form (|14f> has been 
created. This is a great result - it means that we do not need to know the number of atoms in the initial state. So 
long as the total number of atoms is a multiple of three, we get l|14f> . 

We would now like to consider what happens when N ^ 3n. To achieve this, it is helpful to define a measure of 
'cattiness', C, i.e a measure of how similar our final state is to (|14J) . For this, we choose [lflj . 

C = i[P a {N)P fj {N)P 1 {N)] 1/ \ (16) 

This can take values ranging from zero, when the state is very different from (|14f> . up to one, when the state has the 
same form as (1141) . In Figure 6, we have plotted the value of C for the final state (created by the procedure outlined 
above) for a range of values of N. We see that we get a state of the form of 1)140 if N is a multiple of three, but not 
otherwise 20]. This, however, should not overly concern us. For one, this is a very simple procedure for creating cat 
states. The straightforward implementation of this scheme far outweighs the fact that it only works for one in three 
trials, which is certainly not a prohibitively low success rate for usefulness in a range of schemes. Secondly, it is easy 
to post-select trials where the cat was successfully created and throw away all the unsuccessful attempts. 
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FIG. 6: Plot of the 'cattiness' measure C defined by 116H as a function of the total number of atoms in the system, N . A cat 
of the form of l|14fl is created (i.e. C = 1) only if N is a multiple of three. 



IV. DISCUSSION 



One of the great advantages of this scheme over a procedure we previously proposed 0] is that, instead of having 
to control the rate of rotation of the lattice to a high degree of precision, we need only control the time that the 
system is allowed to evolve when the potential barriers are high. This should be much easier to achieve in practice. 
In this section, we investigate how sensitive the cat creation process is to the accuracy of this timing. 

As our benchmark we will take a value of C = 0.9. This corresponds (in the case of the state being symmetric in a, 
/?, and 7) to 90% of the atoms being in a state of the form of ifPUl . We would now like to investigate how accurately 
we need to control the time of the nonlinear evolution in order to achieve a value of C > 0.9 for the final state. If we 
take the time of the nonlinear evolution to be t = (1 + S)2ir/(3U), where S is the error in the timing, we would like 
to find the maximum value of S such that C > 0.9 as a function of N . We denote this maximum value as So- 

A plot of these results is shown in Figure 7. The crosses show the points for total numbers of atoms that are a 
multiple of three and the solid curve is a best-fit line. The approximate scaling from this fit is, 



S fa 0.24/AT. 



(17) 



This means that the error that can be tolerated in the timing scales as At ~ I /UN. We note that this result is 
consistent with the calculation in the previous section for three particles. In i|12|) and (|13|) . the most rapidly varying 
terms are proportional to cos(3[/t), which means that the error in t scales as At ~ 1 / (3C7) . 
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Here we see the great advantage of this scheme: the accuracy with which we need to control the timing scales as 
1/N. This is much more favorable than the exponential dependence on N that the rate of rotation of the lattice 
had to be controlled in a previous scheme ^3] ■ This means that the present scheme has a relatively low premium on 
creating cat states with large numbers and so is particularly well-suited to this purpose. This is good news since the 
advantages of cat states in quantum schemes and investigating fundamental issues of physics are more fully exploited 
the larger the cat state. 

The fact that the non-linear evolution needs to be controlled more accurately the larger the cat is will ultimately 
limit the size of cats that can be created by this process. However, this fact could also be turned to an advantage. 
For example, we could use this scheme to measure the value of U accurately - something that is difficult to achieve 
experimentally One way this could be achieved is to carry out the procedure outlined above for different values of t 
and then make measurements on the final state. Whenever, Ut is a multiple of 2ir/3, a cat of the from of (|14|> will 
be created when N is a multiple of three. This means that one-third of the time we would find all the atoms in the 
same angular momentum mode - a highly unlikely event if a cat state of the form of l|14|) had not been created. This, 
then, provides a clear means of calibrating Ut to within 1/N and, since time can be measured extremely accurately, 
would enable a precise determination of U. 




20 40 60 80 100 
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FIG. 7: Plot of 1/do as a function of N, where So is the maximum value of the error in the timing of the nonlinear evolution, 6, 
that allows C > 0.9. The crosses are numerically calculated data points for values of N that are multiples of three. The solid 
curve is a line of best fit intended to find the scaling. The solid curve should not be interpreted as an interpolation between 
points, since the cat creation process only works when TV is a multiple of three. 



V. READ-OUT AND INTERFEROMETRY 

One interesting application of macroscopic superpositions of flow states is that they may allow for quantum-limited 
measurements of angular momentum, i.e ultraprecise gyroscopes [T^. Il4| . To achieve this, we would like to carry out 
a form of interferometry using cat states. A similar idea makes use of 'quantum beam splitters' to accurately measure 
phase shifts in optical and double- well BEC systems 0, ■ 

The basic idea is to create a macroscopic superposition of supcrfluid flows using the procedure outlined above. The 
cat state is then held for some time At and any rotation of the lattice during this time will shift the energy levels of 
the different angular quasi-momenta and so encode phases on the state. The cat-making procedure is then repeated, 
in analogy with the second beam splitter in a Mach-Zchnder interferometer, and the number of atoms in each angular 
momentum mode is measured. This should enable the rotation of the lattice to be determined accurately. 

In the present discussion, we will limit ourselves to the case where the total number of atoms is a multiple of three. 
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In this case, it can be shown that the transformation of the state by the cat creation process is, 

\N,0,0) afh \ / \N,0,0) a/3j 

\O,N,O) a01 — >CM |0,7V,0) Q/37 1 , (18) 
|0,0,7V) Q/37 / V|0,0,JV) 

where U is the unitary transformation, 

U = ^= I 1 e" 42 ^ 3 I ) . ( I !) ) 




It is straightforward to show that U 3 = 1. This means that, since that cat creation process corresponds to the 
operator U, we pick the second cat operation (corresponding to the second 'beam splitter') to be U 2 . This ensures 
that the second operation is the inverse of the first. Experimentally this could be achieved simply by allowing the 
state to evolve twice as long with the non-linear interactions, i.e. t = 4tt/(3U) p3 |. 

We now need to consider how rotations of the lattice shift the energy of a, (3, and 7. For small changes in the 
angular momentum, AL, we can write, 

AE a — AL = £L, (20) 

where £ is given by L/I and I is the moment of inertia of the atoms. We saw earlier that the angular momenta of a, 
/?, and 7 are 0, h, and —Ti respectively. This means that the respective energy shifts of these modes are 0, £Ti, and 
— and the Hamiltonian describing the system at this stage can be modified from (JHJ to give, 

H = h [-2Ja+a + ( J + + (J - Olh] ■ (21) 

The matrix describing the phase shift of each mode due to evolution with this Hamiltonian for time At, e lHAt / h : can 
be written in the basis {\N, 0, 0) Q| g 7 , |0, AT, 0) a pj, |0, 0, N) a p^} as, 

-i3NJAt g 

Q= I e jJV « A * ) , (22) 
t 




where an overall phase has been neglected. 

Before discussing the results, let us first summarize the interferometry scheme. We first create a cat by the process 
described earlier and then hold it for some time, At. It is during this time that any slight changes in the rotation 
of the lattice are detected. Finally, we repeat the cat creation process, but this time for double the duration of the 
nonlinear interaction, and detect which angular momentum modes the atoms of the final state are in. Even though we 
have two cat-creation processes in this scheme, the overall chance of success in any given trial is still one-third. This 
is because either both processes are successful or both are not depending on whether the total number of particles 
in the system is a multiple of three, i.e. the success of the two processes are not independent. Furthermore, we can 
throw away any unsuccessful trials since we will know that it has been unsuccessful when we measure it if all the 
atoms are not found in the same angular momentum mode, a, f3, or 7. This means that the restriction that the total 
number of atoms must be a multiple of three should not overly concern us. 

The final state created by this scheme, if we begin with the ground state of the system (i.e. all the atoms in the 
non-rotating a mode), is IV") final = U 2 QU\N, 0, 0) a ^ 7 . When £At = 0, i.e. there is no rotation of the lattice or the 
time over which the rotation is applied is vanishingly small, we get Q = 1 and the final state is the same as the initial 
state, as it should be. The probabilities for finding all the atoms in a, 0, and 7 respectively are, 

P a = - [1 + 4:Cos 2 (N£_At) +Acos(N£_At) cos{3NJAt)] (23) 
9 

Pp = - [l + 4cos 2 (iV£Ai- 2tt/3) + 4cos(iV£At- 2tt/3) cos(37VJAt)] (24) 
9 

P 7 = - [1 + 4cos 2 {N£_At + 2tt/3) + 4 cos{N£_At + 2tt/3) cos(37VJAi)] , (25) 
9 

and P a + Pfs + P 1 = 1, as it should when N is a multiple of three. 

We see that the frequency of oscillation of the probabilities (|23|l - (125(1 scales as N. So by changing £At = LAt/I 
by an amount of the order of 1/N all the atoms can be transferred from one angular momentum mode to another. 
This is a dramatic observable that should be able to be measured experimentally, thereby enabling measurements 
of angular momentum with a precision that scales as 1/N. This is the same scaling as the fundamental quantum 
limit (Heisenberg limit) and, since N can be large for macroscopic superpositions, opens up the exciting possibility 
of achieving ultra-precise gyroscopes. 
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VI. CONCLUSION 



We have presented a detailed scheme for creating macroscopic superpositions of different superfluid flow states in an 
optical lattice. Such states are of fundamental interest in exploring the boundary between the regimes where quantum 
and classical physics are valid. The scheme we have proposed has some pleasing features. BECs in optical lattices 
present a system free from impurities that is highly controllable - the coupling between sites and the interactions 
between atoms can be adjusted over several orders of magnitude. Furthermore, this scheme is very simple and well 
within reach of present experiments. Macroscopic superposition states can be created simply by rapidly switching the 
intensity of the laser beams that form the optical lattice. This scheme has the great advantage that there is a relatively 
low premium on creating superpositions with large numbers of atoms. The timing of the nonlinear interaction needs 
to be controlled to an accuracy that scales as 1/N, which is much more favorable than the degree of control required 
in other schemes. This suggests that this scheme may present a viable route to achieving macroscopic superpositions 
of superfluid flow states in the laboratory. Along with enabling detailed studies of quantum mechanics, these states 
may have interesting technological applications including ultra-precise quantum-limited gyroscopes. 
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